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ABSTRACT 


A method for nonlinear system design, based on the generalized 
Popov stability criterion, is developed. It is shown that, although the 
design method is applicable for all values of q in the Popov theory, itis 
particularly useful when q iS non-zero, since it greatly simplifies the 
design effort for this case. All designs are accomplished in conjunction 
with the Nyquist and modified Nyquist loci. Basic to the design proce- 
dure is the development and utilization of the modified frequency-response 
polar loci for lag and lead compensation networks. Three examples, one 


with digital simulation, are included to illustrate the procedure described. 
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[<x INTRODUCTION 

The now well-known Popov criterion has provided a new approach 
to the idea of absolute stability in the dynamics of automatic control 
systems. As aresult, the large amount of subsequent research has 
contributed much to the development of modern stability theory. 

The importance of stability theory is quite evident when one realizes 
that all the current conventional design techniques for linear, time-invar- 
iant systems are directly or indirectly derived from stability theory. This 
is reflected in the fact that the basic engineering design approach is to 
determine how an unstable system can be altered to make it stable. 

When the system is nonlinear or time-varying, the design problem 
becomes much more complex. However, considerable progress has been 
made in this area in recent years. 

This paper considers systems of the standard form of Figure 1. The 
linear element may contain time delays or distributed parameters. The 
nonlinear element may be time-varying and may contain hysteresis. The 
only restriction on the nonlinear element is that its input-output character 
istic must lie within a sector bounded by two straight lines passing 
through the origin. There is no loss of generality if one line corresponds 


to the x-axis since this can be accomplished through a simple transfor- 


mation. Hence: 
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This class of systems may be treated, in some cases, by phase- 
plane methods or by approximate methods such as the single or dual- 
input describing function. Some exact methods are also historically 
available. 

The solution of the absolute stability problem was first formulated 
in 1944 by Lur‘'e and Postnikov [1] for the case when u =f(e) is located 
completely in the first and third quadrants (k =). Also, in their paper, 
the utilization of a Lyapunov function of the type "a quadratic form plus 
an integral involving the nonlinearity" was first proposed. In further 
works, Lur‘e developed a method which permits one to write down a 
finite system of quadratic equations directly from the given differential 
systems. He showed that if this system of equations, which he called 
resolving equations, has real roots, then this is a sufficient condition 
for the absolute stability of the system of Figure 1. 

Yakubovitch [1] further developed this approach by showing that the 
system of resolving equations can be reduced to a single algebraic 
equation. He then proved that the existence of two real roots for this 
equation is a sufficient condition for absolute stability. 

Most of the subsequent research was centered around the application 
of Lyapunov's method for determining system stability. However, up 
until this time, all the exact methods suffered from the disadvantage 
that even the solution of the simple problems became so formidable 
computationally that the methods were generally unusable as engineering 


bools < 


Then, in 1959, Popov proposed an entirely different approach to the 
problem of absolute stability. He expressed his sufficient conditions 
in terms of the frequency response of the linear portion of the system, 
thus giving his criterion a very simple form that is convenient for graph- 
ical application. In fact, the Popov criterion is So powerful that all the 
results described above and connected with the Lyapunov function, con- 
sisting of a quadratic form plus an integral of the nonlinearity, are 
included in the criterion. 

In the Popov method only the Nyquist diagram of the linear element 
and the inequality constraints of the nonlinear element are needed. 
Hence the criterion is readily applied to higher-order systems and those 
containing time delays or distributed parameters. Since itis basic to 
the design method developed in this paper a brief development of the 
Popov theorem follows. Fora proof of the theorem, see [1]. 

In the most simple case, the problem of absolute stability of dynamic 


processes consists of finding when the equilibrium state of the general 


system 
x= Ax + bu 
=e) (2) 
e = C'X 


where x is a vector, AiS a square matrix, bisS a column vector, c' isa 
row vector, and u is a scalar, is on the whole asymptotically stable in 
the Lyapunov sense. The formulation (2) is a system of differential 


equations differing from a linear homogeneous set with constant 
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coefficients only by a nonlinear function whose argument is any one ora 
linear combination of the state variables. 

A, b, and c are real, time-invariant matrices in which some of the 
components may be zero. The nonlinear characteristic, f(e) , is an 
arbitrary, single-valued, piecewise continuous real function, defined 
for all real values of e and satisfying the conditions 


£(0) = 0 (3) 


Gs 10) =e 
e 
which are interpreted geometrically in Figure 2. 

The problem of establishing the absolute stability of the system (2) 
requires the following definitions. A system (2) such that all the roots 
of its linear portion characteristic equation lie in the left half-plane will 
be the so-called principal case of (2). However, a case when some of 
the roots of the characteristic equation are on the imaginary axis and the 
rest are in the left half-plane will be called a particular case of (2). In 
(2), assume that f(e) = he with 0 < h < k , for the principal case and 
with 0 < h < k for the particular cases. Then (2) becomes the linear 
system 

x = Ax + hbe 

(4) 
c=c's, 
Now, Popov expressed his results in the frequency domain. In addition 
to the usual open-loop frequency response, G(jw) (see Figure 1), he used 


a modified frequency response, G*(jw), which is defined as follows 


Il 


u=ke-, 


u = f(e) 


Fig. 2. The input-output characteristic of the nonlinear 
element, 


TZ 


Re G*(jw) = Re G(jw) = X 
(5) 

Im G*¥(jw) = WIm Gijw) =Y. 

From linear theory, for the system (4), obtained from the principal 
case of (2), to be stable for any h inthe interval 0s h <k, itis 
necessary and sufficient that the locus of G(jw), and hence also the 
locus of G*(juw), never intersect the real axis from negative infinity up 
to and including the point (-1/k). Also for the system (4), corresponding 
to particular cases of (2), to be stable with O<h < k , itis necessary 
and sufficient that the locus of G(jw) (or of G*(jw) ) never intersect the 
forbidden zone (the interval of the negative real axis from negative 
infinity to - 1/k) and that the system (4) be stable-in-the-limit, that is, 
that the system (4) be stable for sufficiently smallh > 0 (thatis, 

h = €>0 where ¢€ is arbitrarily small). 

Consider now the original problem of establishing the absolute 
Stability of the system (2), thatis, when f(e) is an arbitrary (not 
necessarily linear) characteristic contained in the sector (0,k) for the 
principal case, or in the sector (€ , k) for the particular cases. 
Additional conditions must be imposed on the frequency response in 
order to guarantee stability for the system (2), not only for linear charac- 
teristics f(e) =he , but for any nonlinear characteristics f(e) contained 
in the given sector. Popov stated these additional requirements in the 
following famous theorem [1]. 

Theorem, For the system (2) to be absolutely stable in the sector 


(0,k) for the principal case, and in the sector (e€ , k) for the particular 


KS 


cases (where e>0 is an arbitrarily small number), it is sufficient that 
there exist a finite real number q such that for all w 2 0 the following 
inequality is satisfied 

Re(1 + jqw)G(jw) + 1/k>0 (6) 
and, additionally for the particular cases, that the conditions for 
stability-in-the-limit be satisfied. 

In order to give the criterion a geometric interpretation let 

Cru) = xy (7) 
Then 

Re(1 + jqw)G(jw) = ReG(jw) - qwIlmG(jw) = X-qY. 

Hence (6) can be written as 

X-qY+1/k >0 forallw20. (9) 
But the equation 

X-qY+1/k=0 (10) 
is the equation of a straight line (the Popov line) with slope 1/q which 
passes through the point -1/k on the real axis. Hence it is possible to 
give a geometric interpretation to the Popov theorem. 

For the system (2) to be absolutely stable in the sector (0,k) for the 
principal case, or in the sector (€ , k) for the particular cases, itis 
sufficient that there exists in the G*(jw) plane a straight line, passing 
through the point -1/k on the real axis, such that the modified frequency 
response G*(jw) lies strictly to the right of it. Also, for the particular 
cases, the conditions for stability-in-the-limit must be satisfied. Some 


examples are shown in Figure (3). 
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Attempts to extend the Popov theorem have met with considerable 
success. Hsu and Meyer [2] have presented a generalized Popov theorem 
which includes the recent extensions. 

Theorem. Consider the basic system of Figure 1 with the linear 
element being output stable. For the system to be both absolutely 
control and output asymptotic for (u/e) e (0/k), it is sufficient that there 
exists a real number q such that for all real w20 and an arbitrarily 
small 6 > 0, the following inequality is satisfied: 

Re ( (1 + jwa)GUw) ) + 1/k 26 > 0 (11) 

The restrictions on qand k depending on the nature of the nonlinear 
element, are: 


1. for a single-valued, time-invariant element: 
if O0<k < ©, then - »< g < & 


ifk=o, thenQ<q <» 


2. for a nonlinearity with passive hysteresis: 
O<k < © and= oq = 0 

3. for a nonlinearity with active hysteresis: 
jes ok =o andevs= qi. oe 


4. for a general nonlinearity (time-varying, and possibly with 
hysteresis): 


O<k <s® andgq=0 
If k = ~in cases (1), (2), or (4) above, then an additional requirement 


must be met. For all t and for every value e_, there must exist a finite 


ay 
Value ue) SUCK Etna 
[Gey =e oma le(t)| < e, (12) 
It can be noted immediately that the term "absolutely stable" has 


been replaced by "absolutely control and output asymptotic". This is 
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mathematically the least restrictive practical stability definition. It 
simply means that the control signal or the output signal eventually go 
eO Zero. 

A valuable engineering tool provided by this theorem is the trade-off 
available between the requirements on the linear element and those on 
the nonlinear element. For example, if the nonlinearity is single-valued 
and time invariant, a non-zero q may be chosen, thus providing a greater 
range for G(jw). 

The Popov criterion provides another valuable asset to the engineer. 
A measure of the rate of system damping is easily incorporated in the 
basic theorem. If the linear element of the system of Figure l is output 
stable of degree a, this simply means that, for a>0O , the system 
output, in response to either an initial condition or an impulse, tends to 
zero faster than the function cee For a<0, it means that the output 
multiplied by aor ultimately goes to zero. 

Now it can be said that if the system of Figure 1 is control asympto- 
tic of degree a, there exists a real number @ such that for every set of 


initial conditions 


i ee ane pidit <r. (13) 


Also, itis output asymptotic of degree wif there exists a real number a 


such that for every set of initial conditions 


© 


| ett) Pat<e , (14) 


° 


In order for the Popov inequality to reflect a measure of damping, the 
linear element must be output stable of degree a, and G(jw-a@) must be 
substituted for G(jw). Hence, for the system of Figure 1 to be control 
and output asymptotic of degree a, the following inequality must be 
satisfied: 

Re ( (lt+jqw)G(jw-a) ) + 1/k 26 > 0 (is) 

Methods of system design based on the Popov criterion have been 
developed by several authors. See, for example, [2], [8], and [4]. 
These methods are all relatively straightforward and easily applied for 
the case q = 0, beacuse of the invariance of the Popov line with frequency 
on the Nyquist plot. However, for the case q # 0, the design of lag-lead 
compensating networks involves considerable design effort, much of it 
devoted to trial-and-error computation. 

A great deal of this labor is eliminated in the new design method 
proposed here for the large group of systems in which design requirements 
dictate using the case q #0. In order to retain the unparalleled utility 
of the polar-plot representation, the lag-lead compensating networks 


are designed on the G*(ju:) plane. 
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II. DESIGN METHOD 

A sample design objective is first illustrated in terms of a particular 
system. Then the general design method is developed, 

Assume that the linear element of Figure 1 is output stable and has 
the Nyquist plot shown in Figure 4. Suppose also that the system design 
requires the nonlinear element to be a saturating amplifier with the input- 
output characteristic of Figure 5. 

The generalized Popov criterion for a general nonlinearity (q=0) gives 
a maximum value of k = 1.67 for the upper bound on the Popov sector to 
guarantee control and output asymptoticity of the system. This is shown 
in Figure 4, However, since the particular nonlinearity of this system 
falls under Case 1, a non-zero q may be chosen to obtain a less con- 
servative result. In order to do this, the simple transformation is made 
from the G(jw) plane to the G*(jw) plane by multiplying the ordinates of 
the G(ju) curve by their respective frequency values. The result is the 
modified frequency response, G*(jw). Assume that G*(ju) has the form 
shown in Figure 6. Drawing a tangent line to the G*(jw) curve (see 
Figure 6) at the pees eae of the G*(je') curve and the negative real 
axis gives a value for q of 0.1 anda value for k of 2.0. 

The system nonlinearity (see Figure 5), however, requires a Popov 
sector with an upper bound of at least 9.0 in order to guarantee a stable 
system with the chosen amplifier. Hence the system must be compen- 
Sated. Excessive compensation and a too conservative solution can be 


avoided by choosing a non-zero qd, The design of lag or lead cascade 
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compensation of the G* plane for systems of this type can be accomplished 
by the method described below with an appreciable decrease in labor over 
previous methods. The compensator is to be inserted in the output of the 
nonlinear element in the system of Figure 1. Before proceeding with the 
compensator design, several preliminary relationships which are instru- 
mental to the design method are given. 

One of the basic difficulties of designing in the G*(jw) domain is the 
fact that 

IGGw)G.Gw) ] * #G*(jw)G .* Fw) (16) 
where G,(jw) is the transfer function of the compensator. In order to 
show later how this difficulty may be minimized, the exact relationship 
is derived. 

G(ijiw) = ReG(jw) + jImG (jw) 


G (jw) = ReG (jw) + jImG,Gw) 


G* (jw) = ReGGw) + jwimG Gu) a 
G_*(jw) = ReG.Giw) + julmG, (jw) 
Multiplying GGw)G, Gw), 
G(ijiw)G. (iw) = ReGGw)ReG. (iw) ~ IMGGw)ImG,( w) 
+ j{(ReG(Gjiw)ImG,(jw) + ReG.Giw)ImG (jw) J = 
“ (GGu)G. iw) ] * = ReGGw)ReG (iw) - ImGGw)ImG.(jw) 
+ jas[ReG(jw)ImG (jw) + ReG ,Gw)ImG(jw)). s 
Multiplying G*(jw)G.*(Gw), 
G*(jw)G.* (jw) = ReGijw)ReG. (jw) - wImG(jw)ImG, (jw) a 


+ jw(ReG(iw)ImG (jw) + ReG.(jw)ImG(jw)]. 


ae 


Comparing (19) and (20), 

IGGw)G (iw) ] * = G*Gw)G.* Gu) + (w®-1l)ImG(ijw)ImG, (iw) . (21) 

The frequency-response plots of the compensation networks are 
essential to the design procedure. The transfer function of a lag compen- 


sation network is given by 


1+ST 
G46) = 


bo ar Where T, = To (272) 


and that of a lead network is given by 


ieee or, 


G(s) = where T, > To. (23) 


lL aote 
The Nyquist plots of the two networks are shown in Figure 7, 
The geometric configuration of the modified frequency responses, 
G.*(jw), is now derived for the lag filter only. The lead filter gives the 


same result only with 7, and tT} reversed. For the lag filter: 


l+w°aT Dy Ga 
~ gt tae Ts 
Ge* jw) = wee ey a) (2.4) 
nl UT 
7 eta 3 = u® (Ts - Ty) 
ee ae a 
nl Ww Ty 


Solving for w* in the X-expression: 
(lL +w* 7," )X = 140° 7, Te 
XS <= 1 ian, 7. 


err,” xX = a” T; Tm =~ 1 - X 


eee 
AX ar Ty T5 


23 


*s1o[d zeTod yIOM\eU UOTJESUSdWOD pest, pue Hey 





oy "Dia 


24 


Substituting into the Y-expression: 


1-xX 
Ty = To 
2 
a — Tai Te : 


aden te) 


Se 


Peete 7) 


Tee TT. Te 
mere - 1,7. + (EX 
Te — 7, Te 


(1 -X) (Tz - T) 
Tees Ty To - TX + Ty 


(1-X) (T2 - T,) 


Ted — Te ) 


(X-1) (T2- aa ) 


T, (T2- Ty) 
Mnenetore; 

Y= 1/1, (X-1) . (27) 
Similarly for the lead filter, 

Y=1/t2 (X-1) , (28) 
which is, of course, the equation of a straight line. This very conven- 
iént result is instrumental to the design method. In order to compute the 
end points of the line (24) is simply evaluated at w= 0 and w=¢,. The 
results for both filters are shown in Figure 8. The frequencies w =1/7, 
and w = 1/Tz occur at the geometric means of the respective polar plots. 

The ideas developed so far can now be applied to the system with 
the modified frequency response shown in Figure 6. Suppose we 


decide to compensate this system with a lag compensator to force the 


a8, 


*SYIOMIOU UOTIESUSdUIOD 
pest, pue bey] 1ojJ sjzojd iejod essuodsei Aouenbsjj peljIpoywys °g ° Hig 


A) 
es 


26 


(MD) x° Day 





TC Gel 1OJ (MI ¥°D 


(1) x “wy 


compensated locus, [G(jw)G.(jw)]*, to intersect the negative real 
axis ata point 2 - 0.1 and at a certain crossover frequency, ww. 
This would give a Popov sector of at least (0, 10) which would guar- 
antee a stable system with the required amplifier, provided that a 
Popov line of non-zero q can be drawn through the point (-0.1,0). 
The procedure to accomplish this alteration of the G*(jw) locus 
is first outlined briefly, then in more detail. The desired negative 
real axis crossover point for the compensated locus is designated as 


P, . Thenthe 1/G*(jw) and the G,*(jw)/P, loci are constructed so 


0 
as to intersect at the point P, where the value of the 1/G*(jw) 

locus corresponds to the chosen crossover frequency, W,, and the 
Pamamectchs, 7, and Ts , are thien computed so that the trequencics of 


the two curves are equal at the intersection point. This construction 


designs the compensator since it guarantees that 


WEA RIG AUP, “at? = Tal, 
or 

G*Vo)G-* We) =P, at w= ow, 
where G*(jw)G.*(jw) and [G.4jw)G(jw) ]* are related by (21). This 
relationship is discussed in detail later. 

To show how easily this is accomplished in practice, the example 
of this section is used. Let P,= -0.1 andw, = 1.0. The first step 
would be to construct the 1/G*(jw) curve. However, this is unnecessary 


since only one point of the 1/G*(jw) curve is actually used, the one 


Ze 


corresponding to uw, = 1.0. The coordinates of this” pointe. P. sameseccume, 
computed and plotted by inverting the complex number read off the G* (jw) 
locus for w, = 1.0. 

The next step is to construct the G,* (jw)/P, plot by simply rotating 
the G,*(ja) plot by 180° and placing the w= 0 point of the compensator 
at 1/P,. The slope of this straight-line plot is adjusted so that the line 
passes through the point P,. Thus P, corresponds to the interseetienmer 
the 1/G*(w) and the Ga(jw)/P, loci, and the frequency at this point 
must be the same, i.€., W= wW,, on both curves. Refer to Figure 9. 

Equation (27) gives the equation of the lag-compensator plot. By 
comparing this with the standard point-slope form of a straight line, it 
can be seen that the slope is 1/T, . Hence 7, can be determined by 
simply computing the slope of the line graphically. Then Tz is 
calculated from either part of (25) using w= u,. The coordinates, X 
and Y, are simply given by the real and imaginary parts of P, /P, 

This is a very easy calculation since the value of the squared parameter 
(in this case, 7, ) is always substituted numerically, whether the 
filter is lag or lead. Using w= ua), in the calculation of T, guarantees 
that the frequency of the G,*(jw)/P, function at P, is the same as that 
of 1/G* (jw) at P, . 

To avoid a sign error when computing P, /P,, it should be noted that 
X and Y are always positive for a lead compensator and X is positive and 
Y is negative for a lag compensator, It should also be noted that there 


are slope limitations on the straight-line compensator plots. Infinite 
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slope corresponds to the pole and Zero which cancel each other while 
zero slope corresponds to infinite pole-zero separation. 

The frequency response curve of interest in the system design is the 
[G(Gja)G. (jw) |] * curve. The plot given by the compensation method 
explained above gives the G*¥(jw)G.*(jw) curve. Equation (21) relates 
these two quantities. Note that since WwW, was chosen to equal 1.0 in the 
above description, the two curves cross the negative real axis at the 
same point. This is not the case in general. 

Assume that the two resulting curves are as shown in Figure 10. 
Here a non-zero q line can be drawn through the point (-0.1,0) and the 
system design is completed. In an actual design problem there may be 
other performance specifications, such as minimum bandwidth, which 
must be met before the design is complete. These additional requirements 
will influence the choice of w, and P, 

In order to generalize the method to any W<, (21) is repeated for 
reference, 

[G(jw)G gw) ] * = G*GuG.* (jw) + (w*~L)ImG(jw)ImG.(jw). (21) 
For simplicity, let (w°-1)ImG(jw)ImG, (jw) =T. It can be noted 


immediately that 


T = 0 when w = 0 and © since ImG,(jw) = 0 (29) 
value for which ImG(jw) = 0. 
For the compensated system to be worsened by the effect of T, (i.e. 


& 


the [G(jw)G.(jw) ] * curve be more negative than the G*(jw)G,* (jw) curve 
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Fig. 10. Modified frequency-response locus for 
compensated system, 
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at the negative real axis crossover), the term T must be negative for 
wW= W.. 
For a lag network, since ImG,(jw) equals zero or a negative value for 
all w, this can occur only under the following two conditions: 
if w. <1landImG(Gw,.) <0 
or if w, >1andImG(Ujw,.) >0. (30) 
For a lead network, since ImG,(jw) equals zero or a positive value for 
all w, T can be negative only under the following conditions: 
if w,. < 1 and ImG(jw,) >0 
orifw.> 1 and ImG(jw,) <0. (31) 
Using the methods outlined above, a rapid design procedure can now 
be stated. Since the subject of compensation is treated in most control 
system design textbooks, no attempt will be made in this paper to sum- 
marize the various helpful hints given by other authors for compensation 
design of linear systems uSing lag and lead networks. Since the design 
method outlined here involves reshaping a linear-system frequency- 
response curve, itis assumed that the designer has a feel for the steady- 
state and transient effects of compensation on linear system performance 
and is cognizant of the advantages and disadvantages of lag and lead 
compensators. See, for example, (5). With this in mind, the design 
steps may be listed as follows: 
1, Plot the polar curves for G(jw) and G*(jw) aS wvaries from 0 
to ©, The G*(jw) plot is easily made by multiplying the 


ordinates of the Gfijw) locus by their respective frequencies. 
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2. From the required design specifications and system compo- 
nents choose the type of compensation (lag or lead) and W,. 

3. Use (30) or (31) and the design specifications for the initial 
cheice ofeP, . 

4. Design the compensator. 

5. Sketch the G*(jw)G,*(jw) locus by evaluating the compensator 
of Step 4 for a minimum number of frequencies. Sketch, if 
necessary, the [G(jw)G, (jw) ] * locus using (21). 

6. Evaluate the design in view of the specified performance 
Criteria, 

If multiple sections are required the process is simply repeated 
using the [G(juw)G (ju) ] * curve in Step 1. Since the design method is 
virtually identical for lead and lag compensators, a detailed discussion 
of lead network design is left to the examples section. It should also be 
noted that the design method introduced above is quite flexible. No 
additional complexity is introduced when a P, is desired off the negative 
real axis. 1/P,is simply designated at the desired point and the 


G.*(jw) line is rotated the required number of degrees, 
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III, EXAMPLES 

To illustrate how the preceding results might be used, the following 
examples are given. 

Example 1. Itis desired to compensate the system shown in Figure 
ll so that the output signal will be asymptotically stable of degree 0.1. 

To simply guarantee output aSymptoticity with the specified type of 
nonlinearity, the minimum Popov sector must be [0,5]. However, for 
the system to be asymptotically stable of degree 0.1, the Popov line 
(for some q) must intersect the negative real axis at a value <-0.2 and 
must not be closer than 0.1 horizontal units to the G*(jw) curve at the 
closest point. 

The G(jw) and G*(jw) curves are shown in Figures 12 and 13. A lag- 
type compensator and anw,= 1.0 are selected. To avoid the possibility 
of a borderline design, a P, of -0.08 is used, as opposed to 0.1. The 
graphical design of the compensator is shown in Figure 14. The 1/P, 
point is at 1/-0.08=-12.5. To locate point P, = 1/G*(jl), note on 
Figure 13 that G*(jl) = -0.4-j0.2 from which 1/G*(jl) = -2+jl as shown. 
The’ parameter values, % and Tz , are computed using the slope"erare 
G.*(jw)/P, plot (see Figure 14) and Equation (25). The resulting compen- 


sator is given by 


eee A Seele 


aol Ss oe 


G.(s) = 


Using a minimum number of points the G*(jw)G.*(jw) locus (see 


Figure 15) is easily sketched. Then, using (21), with (30) as a check, 
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Pig. 12. G(jw) locus for Example 1, 
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Fig. 13. G*(jw) locus for Example 1. 
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Fig. 15. G*(Gjw)G,.*(jw) locus for Example 1, 
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the [G(jw)G.(Gjw) ] * locus can be quickly determined. In (21), the 
(w°-1)ImG(jw)ImG,(jw) term can be evaluated by simply picking the values 
of ImG(jw) off the Nyquist plot of G(jw) and using the expression for Y 

in (25) divided by w for Im G,(jw). 

In this problem, the G*(jw)G,.*(jw) and the [G(jaw)G.(jw)]* loci have 
the same negative real axis crossover at w,= 1 and are virtually identi- 
cal at other points. The Popov line is as shown in Figure 15 guaranteeing 
control and output asymptoticity of degree 0.12. 

Example 2, Consider the nonlinear control system with time delay 
shown in Figure 16. Hsu and Meyer [2] studied this system and obtained 
the following results: 


1. For aclass 2 or class 4 nonlinearity (q = 0), the Popov criter- 
ion is satisfied for 0 <k < 1.85 (see Figure 17). 


2. For aclass 1 or class 3 nonlinearity q # 0), the Popov criter- 
ion is satisfied for 0 < k < 1.98 (see Figure 17). 


Now it is desired to compensate this system to permit a Popov sector of 
[0,5] for a class 1 or class 3 nonlinearity. 

The G* (jw) locus is shown in Figure 17. Lag compensation is chosen 
with an w, = 1.8. Equation (30) quickly shows that if a P, is chosen 
which gives a minimum Popov sector, the initial design will not produce 
an unstable system but will actually give a greater degree of stability. 
Hence a P, of -0.2 is used. Designing the compensator as before (see 


Figure 18) gives the result 
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Fig. 17, G*(jw) locus for Example 2. 
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Therefore: 

1. For aclass 2 or class 4 nonlinearity (q = 0), the Popov criter- 
ion is satisfied for 0 <k < 4.45 (see Figure 19). This example 
illustrates, as a sidelight, how easily the design method may 
also be applied to the case when q = 0, although the advan- 
tages are not quite as great as when gq #0. 


2. Foraclass 1 or class 3 nonlinearity (q # 0), the Popov criter- 
ion is satisfied for 0<k <5 (see Figure 19.) 


The permissible Popov sector will actually be slightly greater than above 
as indicated by (30). If exact sector limitations are desired, these are 
quickly obtained by using (21) and (30) to obtain the [GGjw)G.(jw) ] * 
TOeGuss. 

Example 3. Consider the actuator or indirect control systemmen 
Figure 20. Hsu and Meyer [2] also studied this sytem using the gener- 
alized Popov theorem. Their intent was to determine the critical feed- 
back gain F = F, such that for feedback gains F =F,., the outpiwieswenia! 
y(t) would asymptotically approach a steady-state value dependent on 
the constant input value, r. 

The modified block diagram of this system is shown in Figure 2l. 


This zero-input model is possible because the constant input, Tea 





be combined with the constant-initial-condition-response term with no 
loss of generality. 


Using the modified frequency response function, 


G* (jw) = Ga*Gw) - jF (34) 
where 
| 7 A | 36 - w* 
Ge! ae a eee (35) 
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Fig. 19. G*(jw)G.*(jw) locus for Example 2. 
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For asymptotic stability, the Popov criterion with k = © then requires 
that 

ReG* (j w)-qImG* (ju) = ReGa* (ju')-q (ImMGa* (ju)-F) > 0. 

Go*(ju) is plotted in Figure 22, Since the ImG2*(jw) has a maximum 
value of approximately 0.006, ne Popov criterion is satisfied with some 


q >0 if F > 0.006 and the system is then control and output asymptotic. 


Hence e(~) = 0 and 
G, (0) l 
o) = ee eee 36 
ye) = F460) F+l cr) 


for this system. Therefore the range of F for which the system is 
control and output asymptotic is extremely important since F directly 
affects the steady-state error. 

Digital simulation of the system, with F = 0.006, provided results 
agreeing with the theoretical values. The response of the simulated 
system to an input value of r= 1.0 is shown in Figure 23. 

Assume now that design requirements specify a steady-state error 
which requires the system of Figure 20 to be asymptotically stable for 
F =0.004. Hence the system must be compensated so that the maximum 
value of the ImGs* (iw) is less than 0.004. For comparison purposes, 
the unit step response of the uncompensated system with F = 0.004 is 
shown in Figure 24. Lead compensation, w, = 8.8, and P, = —-(@R7eawe 


chosen. The compensator design is shown in Figure 25 with the result 


leat eves 


Ce) = a guaree (37) 
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Fig. 23. Response of the system of Figure 20 to a unit step 
input. F = 0.006. 
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Fig. 25. Lead compensation design plot 
for Example 3. 
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The G,*(jw)Ge*(jw) locus is shown in Figure 26. 
It should be noted that, in this case, the T term of (29) is of no 
consequence since itis only real and will have no effect on the maximum 


imaginary value of the compensated curve. Hence 


max Im [G.(jw)Galiw) ] * = max Im G,*(juw)G_* (jw) (38) 
and the design is completed. The response of the compensated system 


with F = 0.004 andr = 1.0 is shown in Figure 27. 
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Fig. 27. Response of the compensated system of Example 3 
to a unit stepinput. F = 0.004. 
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IV. CONCLUSION! 

A design method, based on the Popov stability criterion, has been 
developed for the large class of nonlinear systems which can be modeled 
in the configuration of Figure 1. The design procedure uses the modified 
frequency-response polar loci for the linear portion of the system. This 
permits the engineer to utilize the clear physical picture inherent in the 
Nyquist plot. 

As a major step in the development of the design procedure, it was 
shown that the modified frequency-response polar plots for lag and lead 
compensation networks are straight lines. This provided the basis for 
the development of a rapid, straightforward design method which could 


be quite useful in engineering practice. 
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